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Abstract. In this paper we first derive a general method which enables one to 
create expressions for vortices and monopoles. By using this method we construct 
several order-parameters describing the vortices and monopoles of Bose-Einstein 
condensates with hyperfine spin _F = 1 and F = 2. We concentrate on defects 
which are topologically stable in the absence of an external magnetic field. In 
particular we show that in a ferromagnetic condensate there can be a vortex which 
does not produce any superfluid flow. We also point out that the order-parameter 
space of the cyclic phase of F = 2 condensate consists of two disconnected sets. 
Finally we examine the effect of an external magnetic field on the vortices of a 
ferromagnetic F = 1 condensate and discuss the experimental preparation of a 
vortex in this system. 



1. Introduction 

During the last ten years Bose-Einstein condensates (BECs) of alkali atoms have 
turned out to be an excellent system to create and observe several interesting 
phenomena, such as topological defects |2]. The best-known topological defect 
is a vortex, which in a typical single component BEC appears as a long-lived line-like 
singularity in the particle density. In a non-rotating trap, a vortex state cannot be 
energetically the ground state of the system, but its decay is prevented by topological 
reasons. The continuous deformations of the order-parameter which are needed in 
order to reach the ground state require more energy than what is available from 
e.g. thermal excitations. In the presence of dissipation the vortex can move to the 
boundary of the condensate and vanish, but even then it is stable as long as it stays 
in the condensate. The creation of spinor BECs has made it possible to have more 
complicated vortices and other topological structures than what are allowed by a 
single component condensate By spinor condensates we mean BECs which 

have all spin components trapped simultaneously and where spin dynamics between 
different spin components is possible. For these reasons spinor condensates allow 
for richer topological structures than single component condensates. In experiments 
spinor condensates are realized by using an optical trap to trap the condensed atoms. 
If a BEC is in a magnetic trap, only particles which are in a low-field seeking state 
with respect to the quantization axis determined by the local magnetic field remain 
trapped. When an optical trap is used there can still be magnetic fields present. They 
are not needed to trap atoms, but e.g. to diminish the effect of stray magnetic fields 

HIEICIIHI 

The existence of the topological defects is based on the fact that a BEC can be 
described by an order-parameter ip. The order-parameter is a map from some region 
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of physical space into the order-parameter space M. By examining the properties of 
the order-parameter space one can see what kind of, if any, topological defects are 
possible. This examination can be carried out with the help of the homotopy groups 
of the order-parameter space. 

The use of the theory of homotopy groups to characterize the topological defects 
of physical systems was first used during the late fifties [HI, but became widely known 
only in the seventies; see e.g. pHl- Since then it has been successfully appHed in several 
fields of physics, such as condensed matter physics, particle physics and cosmology 
jlOl nil I12L \1?>\ I14L [T^ . Homotopy groups classify maps which can be continuously 
deformed into one another. Homotopy is a mathematical notion giving an exact 
meaning for this kind of deformation. Homotopy groups have turned out to be an 
effective way of characterizing and classifying topological defects. This classification 
can be achieved if the order-parameter space M is identified with a quotient space 
G/H , where G is a group that acts transitively on the order-parameter space and H 
is a suitably chosen subgroup of G. 

If G/H is known, information on the topological defects can be obtained by 
calculating the homotopy groups TTn{G/H), n = 0,1,2,3. From these tto{G/H) 
characterizes domain walls, 'Ki{G/H) vortices and one-dimensional non-singular 
defects, tt2{G/H) monopoles and two-dimensional non-singular defects, and ^^{G/H) 
three-dimensional non-singular defects. Domain walls, vortices and monopoles are 
defects where the order-parameter has to vanish at some point of the physical space. 
The non-singular defects are defects where the order-parameter is nonzero everywhere, 
and the topological stability is imposed by the boundary conditions. The elemets of 
the homotopy group label the order-parameters in such a way that those labelled by 
the same group element can be continuously converted into one another, whereas if 
the configurations are labelled by different group elements, this is not possible. If 
iTn{G/H) = {e}, i.e. the nth homotopy group is a one element group, no topologically 
stable defects characterized by 7r„ (G/iJ) are possible. 

This paper is organized as follows. In section 2 we derive a systematic way to find 
expressions for vortices and monopoles. This method is based on the properties of the 
relative homotopy groups, and it has not been presented before. In section 3 we review 
the properties of spinor condensates and their ground states. As a new result we show 
that the order-parameter space of the cycHc phase of = 2 condensate consists of 
two sets which are disconnected. This is in contrast to the order-parameter spaces 
of other ground state phases, which consist of one connected set. In section 4 the 
method derived in section 2 is applied in the context of spinor BECs to create order- 
parameters describing vortices and monopoles. Most of the expressions for defects have 
not been presented before. These are then used to find the minimum energy states 
of the defects and to study the superfiuid velocity and angular momentum induced 
by them. In particular we show that in a ferromagnetic condensate the presence of 
a vortex does not have to lead to superfiuid fiow. In section 5 we show that the 
vortex of a ferromagnetic F = 1 condensate derived in section 4 can exist also if the 
conservation of magnetization is taken into attention. Additionally we show that in 
a ferromagnetic F = 1 condensate the conservation of magnetization may lead to 
stabilization of defects which are not stable if the magnetization is allowed to vary 
freely. In section 6 we propose a way to create vortices in a ferromagnetic condensate. 
These vortices are stable in the absence of an external magnetic field. This way is based 
on the use of a topological vortex creation method used before to create defects which 
are stable in the presence of an external field. Finally in section 7 we summarize the 
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results of the paper and in appendix A we give the spin and rotation matrices needed 
in the calculation of defects. 

2. Finding expressions for vortices and monopoles 

2.1. Introduction 

In this section we develop a method for finding explicit expressions for vortices and 
monopoles. The content of this section is quite mathematical, and those interested 
only in the applications may jump to the next section. We briefly introduce some of the 
necessary concepts and flx the notation. To save space we do not deflne all mathematial 
conceps needed; an interested reader can find the definitions in any standard book on 
homotopy theory, such as ^6]. We also omit all proofs here; those can be found in 
jlOl I17j . From now on all maps are assumed to be continuous. In the following G 
is a Lie group and is a closed subgroup of G. The subgroup H can be written 
a.s H = Ho U Hi U H2 U ■ ■ ■, where i/^'s are the path components of H. The path 
component containing the identity element e is denoted by Hq . It is a normal subgroup 
of H and thus the quotient space H/Hq is a group. The path components of H are 
cosets of H in Hq, i.e. Hk = hkHo for some hk £ H. We define the n-dimensional 
disc by Z?" = {x e M" | |x| < 1 } and the n-sphere by 5" = {x e R"+i | |x| = 1 }. We 
denote the nth relative homotopy by 7r„(G, H, e) and an element of this group by [/], 
where / : (£)", 5*"^^, sq) ^ {G,H,e) is a map. Here the notation means that D" is 
taken to G, S"^^ to H, and sq G 5'""^ to e by /. By [/] we denote the equivalence 
class determined by /. It consists of all maps which are homotopic via maps of this 
type. If H — e, we write 7r„(G, e) = 7r„(G, e,e). If G is path connected, 7r„(G, 51) 
and 7r„(G, 32) are isomorphic for all gi,g2 € G, and we use the notation 7r„(G) to 
denote any 7r„(G, gi). There is an exact sequence of homomorphisms between relative 
homotopy groups. This sequence reads 

^ 7r„(G, e) ^ 7r„ (G,iJ,e) ^ 7r„_i(i7, e) ^^V„_i(G, e) • • (1) 

I P* 

^n{G/H,H) 

We have included in the sequence , which is an isomorphism determined by the map 
p:G^ G/H,g^ gH. If [/] e ^„(G, H, e) then = [pf] G ^„(G/ff, H). 

2.2. Physical applications 

Next we consider how previous results can be applied in physical systems. Now G is 
a group that acts transitively on the order-parameter space M. The action of g e G 
on X e M is denoted hy g ■ x. We now choose an arbitrary element Xref & M, called 
the reference order-parameter, and define H = {g Q G\g ■ Xref = Xref}- This is 
the isotropy group and it is a closed subgroup of G. The order-parameter space M 
can then be identified with G/H. The correspondence between the elements of M 
and G/H is a; <^ 9H,g ■ Xref = x. If A is a topological space and f : A ^ G is 
a map, then the map pf : A ^ G/H gives a map from A to the order-parameter 
space, which is now represented by G/H. If the order-parameter space is represented 
by M the corresponding map from A to M is given hy c : A M such that 
c(a) — f{a) ■ Xref for all a G A. When discussing the physical applications we assume 
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that TT2{G,e) — TTi{G,e) — T:o{G,e) — {e}, when 71 and 72 become isomorphisms. 
These conditions hold for M and SU{2), which are the groups used in this paper. 

2.2.1. Monopoles Because 7ri(7J, e) = 7ri(_ffo,e), we see that Tr2{G,H,e) and 
'!Ti{HQ,e) are isomorphic via 72. The map 72?^^ gives an isomorphism between 
'!T2{G/H, H) and 7ri(7Jo,e). This isomorphism can be used in the calculation of 
iT2{G/H, H), since usually it is quite easy to see what 7ri(i?o, e) is pi llflllTT] . Assume 
that [/] £ T^2{G, H, e) and that /[^i is the restriction of / to S"^, the boundary of D'^ . 
Then [f\gi\ G 7ri(-ffo,e). The isomorphism 72 is given by the map [/] 1-^ [/jsi]- We 
use {9, ip) as the coordinates of D^, that is, 9 is the distance from the center of the disc 
and ifi is the azimuthal angle. We take the radius of to be tt. Let g : [0, 27r] Hq be 
a map for which g{0) = .g(27r) = e and let g : — > G be such that g{9 = tt, <p) = g(ip). 
If we also define so = 0), then [g] £ tt2{G, Ho,e),\pg] £ 'JT2{G/H, H) and 
[5I51] = [ffle^^] = [.9] •= 7ri(iJo,e). As explained above, [pg] and [g] are mapped 
into each other by isomorphism 72^^^- 

Next we construct a map from physical space into the order-parameter space G/H 
which describes a monopole with some given winding number. Now we use spherical 
coordinates (r, 9, (p) as the coordinates of the physical space M'^ and assume that the 
monopole is located at the origin. The monopole we construct is independent of the 
r-coordinate. This assumption is not necessary, but to avoid further complication we 
use it here. Let >1 = \ {qi._ ^gg^^g f : A ^ G/H by fir,9,ip) = g{9,ip)H. 
Then for each fixed r > [f] £ n2{G/H,H) and [/] is the unique inverse image of 
[g] £ 7ri(7?0)e) in the map 72P^^. These elements have the same winding number. 
Thus if one wants a map f : A G/H, which describes a monopole with a given 
winding number, he has to find g : G such that [ff is an element of 7ri(i?0: e) 

with the wanted winding number. Then / is obtained by defining /(r, 9, ip) = g{9, (p)H. 
The corresponding order-parameter x : A — > A/ is defined by x{r, 9, (p) — g{9, ip) -Xref- 

2.2.2. Vortices The group structure of tto{H, e) has to be defined a little differently 
than that of other relative homotopy groups JOI- We define no{H,e) = H/H^ = 
{Ho,Hi,H2,...} and = [0,27r]. Let g : ^ G be such that 5(0) = e and 
(?(27r) £ Hm for some m £ {0, 1,2,.. .}. If we choose so = then [g] £ 7ri(G, iJ, e). 
The isomorphism 71 : 7ri(G, iJ, e) 7ro(i?, e) is given by [g] i— > [g(27r)], where 
[g{2n)] = Hm. 

We use cylidrical coordinates (r, z,ip) as the coordinates of physical space and 
assume that the vortex is located on the z-axis. We define = \ Me^, where Ee^ 
denotes the ^-axis. We define f : A ^ G such that /(r, z,0) = e and /(r, z, 27r) £ Hm- 
Then for each fixed (r, z) [f] £ 7ri(G, H, e), [pf] £ tti{G/H, H) and the image of [/] in 
the isomorphism 71 is [/(r, z,27r)] = H„i £ nQ{H,e). Thus / gives a vortex with the 
winding number represented by Hm £ T^oiH^e). The corresponding order-parameter 
a; : A ^ M is defined by x{r,z,ip) = f{r,z,ip) ■ Xref- For a vortex constructed this 
way x{r, z, 0) — Xref for all r > 0, z e M. This requirement can be relaxed, but for 
our purposes that is not necessary. 

3. Spinor Bose-Einstein condensates 

In this section we review the ground-state order-parameter spaces of spinor BECs. As 
a new result we show that the order-parameter space of the cyclic phase consists of 
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two disconnected sets. A spinor in one set cannot be continuously converted into a 
spinor in the other set while staying in the order-parameter space all the time. Only 
the order-parameter space of the cyclic phase has this structure, since in other ground 
states the order-parameter space is connected. 

A spinor condensate of atoms with hyperfine spin equal to i^, F = 1, 2, . . ., is 
described in the mean field theory by the order-parameter t/', which can be written in 
the form 7/)(r) = •\/n(r)5(r), where n(r) is the particle density, ^(r) is the transpose of 
the complex vector (^^(r), Cf-i(i"), ■ • • ,C-f(i")) and ^(r)'''^(r) = 1 0^]. In the rest 
of the paper we call ^ the spinor. In some publications the term spinor refers to tp, so 
one must be careful with the terminology. When we determine the order-parameter 
space, we set the density equal to one, so the order-parameter is just the spinor 
If the normalization of ^ is the only restriction imposed on the order-parameter, the 
order-parameter space becomes 5"*^+^. This space allows for no topological defects 
characterized by 7r„ with n = 0, 1, 2, 3. This is because from the theory of homotopy 
groups it is known that 7r„(S'^^+^) = {e} for n ~ 0, 1, 2, 3 and F = 1, 2, 3, .... However, 
for example the order-parameter space of an F = 1 condensate may be only a subset 
of . This can be inferred from the energy functional, which for an F = 1 condensate 
in the absence of external magnetic field reads ^ EI 

= f dh{^ v^*(r).V^.(r) + y(rHr) + ^K+/3i(F)2]}. (2) 

•' i=-l 

Here F is the (hyperfine) spin operator, (F) = ^'l'(r)F^(r), M is the mass of the atom, 
V is the external potential, ai = = and ap is the s-wave 

scattering length in the total spin F channel. 

One sees that it is energetically favoured that either |(F)| = or |(F)| — 1, 
corresponding to the cases /3i > and (3i < 0, respectively. When (3i > the 
system is said to be antiferromagnetic, whereas ii (3i < the system is said to be 
ferromagnetic. It turns out that a group that acts transitively on the set of spinors 
fulfilling the condition |(F)| = 1 is SU{2). This can be shown using a similar reasoning 
as that shown below in the context of the ground-state phases of F = 2 condensates. 
SU (2) acts via its irreducible three-dimensional representation, i.e. the spin rotations 
of an F = 1-particle. For spinors with |(F)| = we choose G = M x SU{2), where M 
gives the gauge transformations of the spinor as ^ i-^ e*^^, 6* € M. As before, SU{2) 
describes spin rotations. 

t In principle the order-parameter space should be written as M' = (R+ x M) U {0}, where 1R+ is the 
set of real numbers larger than zero giving the possible values of the square root of the density, and M 
gives the order-parameter space related to the density-independent part of the order-parameter. Now 
M consists of the possible values of the spinor, while in the case of a single component condensate 
M = , which characterizes the possible values of the phase of the order-parameter. The point 
{0} denotes the case where the density is zero. One sees that if M' is the order-parameter space 
there are no topologically stable defects. Any order-parameter can be converted into any other order- 
parameter via deformations which reduce the density to zero in an appropiate region of the physical 
space. Prom a physical point of view this is unlikely to happen, since reducing density to zero is 
not energetically favorable. In principle the necessary energy could come for example from thermal 
excitations, but in practise this is unlikely to happen. Thus one can ignore the zero of density. Then 
the order-parameter space becomes R+ x M. For this 7r„(R+ x M) = 7r„(R+) x 7r„(M) = 7r„(M), 
since 7r„(R+) = {e}. Thus, from the point of view of topological defects, it is enough to study the 
structure of M only. 
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Here a2 = y(4g2 + 3.94), /32 = -7(52 - 9i), and 72 = ^(go - .94) - 1(32 - gi), where 
^. ^ 47rn^a. Q ^ 2^2^-2 - 2^i^_i + . As in ail F = 1 system, the energy 
is invariant in position-independent spin rotations and gauge transformations. The 
possible ground states have been calculated in and can be classified as follows, 

(i) If /32,72 > the energy is minimized when |(F)| = 8 = 0. Spinors with these 
properties are called cyclic, (ii) When /?2 < 0, 72 > the minimum is obtained by 
making |(F)| = 2, 9 = 0, and the ground state is ferromagnetic, (iii) If /32 > 0,72 < 
the minimum is achieved by maximizing 6, i.e. |9| = 1, and |(F)| = 0. The ground 
state is called polar or antiferromagnetic. (iv) Finally, if a2 and (32 are both negative, 
the ground state is ferromagnetic for 4|/32| > I72I, and polar otherwise. This is because 
1(F) I and |0| cannot be maximized simultaneously. 

Next we discuss briefly the structures of the order-parameter spaces of the 
ground states. Since the interaction energy is invariant in gauge transformations 
and spin rotations and /32,72 are arbitrary, also |9| and (F)^ are invariant in gauge 
transformations and spin rotations. Furthermore, for every spinor ^ there is always a 
spin rotation R which rotates the spinor so that in the rotated state the spin is 
parallel to the z-axis, i.e. (F) = {rIYFR^ = {RQ^fP^R^e^ = {F^)e^. In the rotated 
state (Fx) = {Fy) = 0. Thus the spin-dependent terms in equation (j^J can be written 
as l32{Fzy +72|9p. In the cyclic phase the energy is minimized when (F^) = |9| = 0. 
Numerical calculations show that the solutions to these equations are (up to a 
gauge transformation and a rotation about the z-axis) CO = ^(1,0,-^/2,0,-1)"^, 



CI = ^(1,0,0, V2,0)^ and CI' = ^(0,^2,0,0,1)^ (211 121 • All cyclic spinors 



can be obtained from these spinors by a gauge transformation and a spin rotation. 
By exploiting the rotation matrix shown in the appendix one can show that CI, CI' 
can be rotated into each other, while for CO, CI and CO, CI' this is not possible. This 
means that the order-parameter space of the cyclic phase consists of two disconnected 
sets, a fact that has not been pointed out before. Formally this can is expressed as 
■kq{G/H) = Z2. The group G = K x SU{2) acts transitively on both disconnected 
sets. In the ground-state the system can consist of regions, some of which are in a 
CO state and others are in a CI state. These regions are separated by a domain wall. 
This kind of structure is not possible in other zero fleld ground-state phases of F = 1 
or F = 2 condensates. 

In the case of ferromagnetic ground states, the situation is simpler. Now the 
equations to be solved are |9| = {Fx) = {Fy) = and \ {Fz) \ — 2. The only solutions 
(up to a phase) are |F = 2, = 2) and |F = 2,m_F = —2). These spinors can be 
rotated into each other, so the order-parameter space is now connected. By examining 
the spinors obtained by a spin rotation from the reference order-parameters one sees 
that in the case of a ferromagnetic condensate we can choose G = SU{2) instead of 
R X SU{2). 

Similar study for the polar phase shows that the order-parameter space of the 
polar phase is connected but larger than R x SU(2). This complicates the study of 
polar defects, and they will not be discussed here 0. 

Before going to the details of defects we review the possible ground-state phases of 
some alkaH atom condensates. Experimental and theoretical results indicate that the 



For an F = 2 condensate the energy functional is p3 
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F = 1 spinor Bose-Einstein condensate of ^''Rb is ferromagnetic ^EIIHIEHI- On the 
other hand, the F = 2 condensate of ^"^Rb is probably polar jElElEBI. The ^^Rb F = 2 
condensate seems to be polar jJH] and the ^^Rb isotope with F = 2 ferromagnetic flSj. 

^■^Na scattering lengths determined in indicate that ^^Na F = 1 is 

antiferromagnetic, as has been predicted in ^ and seen in experiments [^. The ground 
state phase of ^^Na F = 2 spinor condensate appears also to be antiferromagnetic 
Experimental study of this condensate is difficult because the |F = 2,m_F = 0) state 
decays within milliseconds pBJ.. 

4. Examples of defects in spinor condensates 

Next we present several examples of vortices and monopoles in spinor condensates. In 
these examples cyHndrical coordinates (r, and spherical coordinates {r,6,(p) are 
used when discussing vortices and monopoles, respectively. In vortices the coordinate 
dependence is such that at <p = the reference spinor is obtained. The vortex cores are 
assumed to be straight and located on the z-axis. Monopoles are located at the origin 
of the spherical coordinates. In what follows the minimum energy of a vortex with 
given winding number is also studied. This is calculated under the assumptions that 
the vortex stays fixed at the z-axis and the density n is independent of (p. The ground 
state density can be assumed to be cylindrically symmetric if the external potential 
used to confine the condensate is cylindrically symmetric and the vortex does not 
split. The latter is not true in general, since usually if a condensate is trapped in a 
harmonic trap it is energetically favorable for a vortex with winding number larger 
than one to split into winding number one vortices. In this case the particle density 
cannot be cylindrically symmetric. The splitting can possibly be prevented by using 
a trap which is steeper than harmonic jJE] or by applying a repulsive potential in the 
vicinity of the rotation axis j22], p 245. 

4-1. Ferromagnetic condensates 

Next we briefly examine the vortex conflgurations of ferromagnetic condensates. 
Here we mean by ferromagnetic condensates systems whose ground state spinor is 
\F,mF = F) or any other spinor obtained from this by a spin rotation. A separate 
gauge transformation is now unnecessary, as a spin rotation alone is able to produce 
that. We choose \F,mp = F) as our reference order-parameter. The isotopy group 
is isomorphic with Z2F and the order-parameter space is SU{2)/Z2f', see j2Sj. The 
isotropy group for F = 1 and F = 2 has been calculated explicitly in [S| . The winding 
numbers of topologically stable vortices range from 1 to 2F — 1. If the winding 
number of a vortex is m, the winding number of the antivortex is 2F — to. Especially, 
a vortex with winding number F is also its antivortex. If F = 1, the order-parameter 
space is SO{3), since SU{2)/Z2 = 5*0(3) [4j. Because 7r2(5f7(2)/Z2F) = {e}, in 
ferromagnetic condensates monopoles are not topologically stable. Next we construct 
order-parameters describing vortices in ferromagnetic F — 1 and F — 2 condensates. 
The rotation matrices needed in the calculations can be found in appendix A. 

4.1.1. Ferromagnetic F — 1 condensate Now H = {I, —1} e SU{2) and, as explained 
in section 2, a vortex is given by a map / : R^^ \ Re^ — > SU (2) such that /(r, z,0) =1 
and /(r, z, 2n) = —I for every 2; G R, r > 0. From the ^/-matrix IIA.4|l one sees that this 
kind of map is obtained by choosing r such that T(r, z,0) — and T(r,z, 27r) — 27r. 
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The functions a and (3 can be arbitrary functions of position. The corresponding 
order-parameter is i/'(r) — \/n(r)^(r), where i^(r) = /(r) • ^^e/ = V'''^\v)£^re f and 
is given in equation ljA.6|l . We get 

((cos ^ — zoos/? sin \ 
- V2e'" sin (3 sin § (i cos § + cos /3 sin §) ■ (4) 
-e^^" sin^ § sin^ /? / 

From the way this was derived it follows that this is defined only for r > 0. This order- 
parameter can, however, be extended to whole by requiring that n(0, z, ip) — 0. 
This is needed in order to keep the order-parameter well defined at r = 0. 

Next we try to find a spinor giving a vortex with the minimum energy. Now 
the superfluid velocity is v = —i^S^^'WS^. When the vortex has the smallest possible 
energy the velocity in radial and z-directions can be assumed to vanish and thus the 
spinor is a function of ip only. Additionally, the density is taken to be cylindrically 
symmetric. Using these assumptions the Euler-Lagrange equations obtained from the 
energy functional show that the vortex energy is minimized when a = constant, /3 = f 
and r(r, z, ip) = ip. Thus the vortex takes the form 

((cos ^ — i cos/3 sin \ 
- V2e*" sin 13 sin f {i cos | -t- cos /3 sin f ) • (5) 
-e^*" sin^f sin^Z? / 

Here we have not set /? = ^, but have allowed it to have any value. This recalls a 
little the expression for the coreless vortex shown in These are, however, different 
defects, since the vortex of equation J^J cannot be converted into a non- vortex state by 
continuous deformations, whereas for the coreless vortex this can be done§. Another 
difference is that in the case of a coreless vortex the particle density does not vanish, 
while in the case of JSJ the density vanishes on the z-axis. The superfluid velocity for 
the vortex state Q is v = —jj^^^^ and if (3 = ^ it vanishes. Thus, in contrast to a 
single-component condensate, in a spinor condensate the existence of a vortex does not 
have to lead to a non-zero superfluid velocity. The same phenomenon can be seen also 
in the orbital angular momentum. For the order-parameter it is L = —Nhcos f3 Bz, 
which also vanishes if (3 — f . The condensate does not have to contain any angular 
momentum although there is a vortex in it. If the condensate is in a state with /3 = 0, 
it contains one unit of angular momentum per particle and the system is similar to a 
single-component condensate with a vortex. In the presence of dissipation the angular 
momentum does not have to be conserved, and the vortex can evolve towards the 
ground state where the angular momentum vanishes. The kinetic energy of © is 
proportional to 3 -I- cos 2/3, which decreases monotonically as (3 increases from zero to 
^ . Thus there is no energetic barrier which could render the vortex state with /3 = 
metastable against conversion into the ground state vortex. 

A vortex obtained from ^ by setting (3 = has been presented before in |4|. 
The general expression for the vortex shown in equation (HJ or the behaviour of the 
angular momentum and superfluid velocity have, however, not been discussed before. 
This also holds true of the ferromagnetic F = 2 vortices discussed next. 

4.1-2. Ferromagnetic F = 2 condensate The isotropy group is H = 
{I, —ioz, (— icr^)^, {—iazY}, where Uz is the z-component of the Pauli matrices. A map 

§ The conservation of magnetization may change the situation. This is discussed in section 5. 
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fm ■ R'^ \ SU{2) represents a vortex with winding number m if /m(r, z, 0) = I 

and fm{r, z, 2tt) — [~iaz)™' for every z g R, r > 0. If we parametrize SU{2) matrices 
as in IIA.411 . these conditions are fulfilled if T(r, z,0) = and T(r, z,27r) = rmr. 
Additionally, for m = 1,3 the condition /3(r, z,27r) — 2fc7r has to hold. Here k is 
an arbitrary integer. For m = 2 the function (3 can be arbitrary. The order-parameter 
becomes 

(cos ^ — i cos/3sin 0^ 
2e''"sin^ sin/3(i cos ^ + cos/3 sin 
^e^*" sin^ (i (cos (3 — cos t cos /3 + i sin r)^ 



tjj{r,z,(p) = ^yn{r,z,(p) 



-3 ^ sin^ 



2 ^-lii (3 {i cos ^ + sin ^ cos /3) 



(6) 



\ e^*" sin"* § sin"* /3 / 

The minimum energy of a vortex with winding number m is obtained when T{r, z,ip) — 
^ and £, = ^ {(p) . Ifm — 2 the vortex minimizing the energy is obtained from equation 
101 by replacing t with ip and setting /S — ^. In this state, the superfluid velocity 
and angular momentum vanish. We show that this is a general property of a winding 
number F vortex in a ferromagnetic condensate with hyperfine spin F. This vortex 
is represented by the element —leH and a vortex can be written as tpi^: v) = 
\Jn(r, z)^(iy9), where t,{'p) = e""'' "'^^j-e/, n = (cosasin/?, sinasin/3, cos/3) is a 
constant vector and ^re/ = \FTmp = F). If /5 = -I we get n = (cosa,sina,0) 
and v^,L, ~ ^t^e = -i^re/e^"^"-^!! • Fe"*'^ "-^ere/ = -lir.!^ ■ F^ref - 0. 

If m = 1, 3 the function (3 that minimizes the energy of a vortex is more difficult to 
find. The obvious choice (3 = is not the correct one, since for example l3{ip) = §1±£ 
produces a smaller energy. Finding the ground states of vortices with m = 1,3 will 
be left for future publications, and will not be discussed here in more detail. If /3 = 
the vortex is simply ip{r, z, ip) = yjn{r, z)e™'^|F = 2, mp = 2), m = 1, 3. 

4-2. Antiferromagnetic F = \ condensate 

It is advantageous to use the U matrix of equation IIA.311 and corresponding 
representation matrices when discussing the defects of an antiferromagnetic _F = 1 
condensate. The reference order-parameter is chosen to be \F — l,mp = 0) and the 
isotropy group is H = {(m2Tr, [/(v?, 0,0)), ((m + ^)2tt, gU{ip, 0,0)) \ (p e [0,47r],m e 
Z}, where we have defined g — U{0, tt, 0). The connected component of the identity is 
Hq = {(O, U{(p, 0, tp)) \ (p £ [0, 2tt]} 3j. Now vortices and monopoles are possible and 
both of them are classified by integers. A general antiferromagnetic spinor can then 
be written as ip{r) = ^n(r)^(r), where ^(r) = e'^*^'''?/*^^) (r)^re/ and C/^^^ is shown in 
l|A.5|l . ExpHcitly 

/ _e-*"W-ijsin/3(r) \ 
V'(r) = y^e'^^"^) cos/3(r) . (7) 

^ e^"W-^sin/3(r) J 

4-2.1. Vortices As can be seen from the isotropy group H, vortices with winding 
numbers m and m -\- ^ are possible (m G Z). Using an analysis similar to that of 
ferromagnetic condensates, a vortex with a winding number m can be shown to be 
given by 10 if 6* fulfills the condition 6'(r, z,0) = 0, 0(r, z,27r) = 27rm. Other angles 
can be chosen freely as long as ^(r, z,0) = ^(r, z,27r) — \F — l,mp = 0). For 
condensates the requirements are 6'(r, z, 0) = 0, 6'(r, z, 27r) = 27r(m + i) and 
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^(r, z,0) — — $(r, z,27r) = \F = l.mp = 0). The superfluid velocity is v = w^^^v 
and it cannot vanish everywhere if there is a vortex in the system. The orbital angular 
momentum is L = Nhl , where I is either m or m + i . In the following we assume 
that Tp{r,z,ip) = ^ n{r, z)£^{ip). The spinor minimizing the energy of a vortex with 
winding number m e Z is obtained when a, (3 are constants and 9{ifi) — mip. If we 
choose = it becomes ipmir, z, ip) — y^n(r, z)e™'^|F = l,mp — 0). A winding 
number m, + i vortex with minimum energy is given by 10 under the conditions that 
0(r, z, 27r) = (2m + i)7r, a is constant and /3 = ^, 

/ -e-^"^sinf \ 
^„+i(r,z,^) = V^^e^("+^)W cosf . (8) 

Properties of vortices of the type l|Hl have been discussed before for example in 
|29l 13(11 13T1I32]. 

4-2.2. Monopoles As explained in section 2, an expression for a winding number 
m monopole can be found as an extension of a map gm ■ [0, 27r] — + Hq, where 
gm determines the element of 7ri(7?o,e) with winding number m. Now we define 
gmiv) = iO,U{amip),0,am{ip))), where am is such that Q;m(0) = 0, am(27r) = 27rTO. 
We define a map gm '■ ^ M x SU{2) such that 5m(7r, = gm{p)- Then fm ■ \ 
{0} [RxSU{2)]/H, fm{r,0,ip) = gm(^, V')^ gives a monopole with winding number 
m. A general form for is gm(0,(p) = {Sm{d,(p),U{am{0,(p), l3m{d,(p),am{0,(p))) 
with the conditions amiT^jf) = ctmiv), /3m(7r, = (5„i(7r, (p) — 0, which follow from 
gm{0 = TT^ip) = gm{^)- From the continuity of it also follows that 5m{^-,p) ~ 
and /3m (0, (p) = (2k + 1)tt with k an arbitrary integer. 

By a simple change of basis equation 10 can be cast in the form 

(cos am sin/3fe \ 
sinam sin/3fc . (9) 
- cos Pk ) 

Here we have substituted 5,am,l3k for 9,a,P, respectively. If (5 = the above spinor 
is a real and normaHzed three-component vector. Then there is an integral equation 



4Wo Jo '^''^■[to^'di) 
which gives the winding number w of the monopole; see ^2] or Using this 

equation one can confirm that by assuming Sm = and that gm fulfills the conditions 
stated above, the winding number of the spinor in equation ^ is m. For a non- 
constant 6m the winding number is the same, as it does not depend on the form of 
6m- From equation one sees that the superfluid velocity vanishes if and only if 6m 
is a constant function. This follows from the fact that 6 = constant <J=4> f^Vf G M^. 
Because v = — i-p-^'''V^ € M^, it follows that v = if and only if 6m is a constant. 

While it is instructive to know what are the general requirements for a spinor to 
represent a monopole, it is important to see that the typical expression for a monopole 
can also be obtained as a special case of above equations. If we choose gm such that 
6m = 0, am{0, f) = niLp and /3m(6', <-p) = {2k + \){-K — 9), where k is an arbitrary integer, 
all the requirements for 6m, ctm and /3m are fulfilled and the monopole becomes 

cos(m(^) sin[(2fc + 1)6* 
^Pm{r, 6) = ^n{r, 9, ^) \ sin(TO^) sin[(2fc + 1)9] | . (11) 

cos[(2fc + 1)9] 
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Choosing fc = gives the usual expression for a monopole. It has been presented 
in the context of BEC before for example in jSH EH, but the requirements for a 
general expression of a monopole and its properties have not been discussed before. 
A monopole with m = 1, fc = has been numerically studied in [Ml j . 



4-3. Cyclic states 

The defects of the CO phase have already been discussed in f3^, so here we concentrate 
on the defects of the CI phase. The isotropy group H has not been calculated before, 
so we have to do it here. For the S'C/(2)-matrices we now use the {/-representation 
shown in ljA.3|l . The reference spinor is chosen to be CI = -^(1, 0, 0, -y/S, 0)^ and a 
general expression for a spinor is 



^(r,z,(p) = ^yn{r,z,(p)e 



id 



V6 



V3 



cos 



4 13 



,2 §. 
1 



[e'^''^(cos/3- cos(2/3)) 



\/2e3*T sm- 

y2cos2 f sin/3] 



sin /3] 

2^, 



V 



^sin/3[-4e3»T cos/3- 
-[e3*''(cos/3 + cos(2/3)) - 

— 



^ \/2sin/3] 
V2 



sm 



2 a 



sin (i\ 



%/3 



-[sin'' I + V2e3*'' cos2 | sin/3] 



.(12) 



Equating this with the reference spinor CI shows that the isotropy group is H = 
{(^(2n + 3m),?7(^,0,0)) | m € Z, n = 0, . . . 5} C M x SU{2) and Hq = (0,1). 
This means that only vortices, not monopoles, are topologically stable. H/Hq = 
H is isomorphic to the group Z x Zg, the isomorphism is given by the map 
(^(2n + 3TO), C/(^,0,0)) 1-^ (m,n) e Z X Zg. Therefore vortices are classified 
by two winding numbers (m, n) and those with different m cannot be converted into 
one another. However, for vortices with winding numbers {m,n) and (m,rt + 6) this 
is possible. 

A vortex with winding numbers {m,n) is given by fm.n ■ \ Ke^ ^ M x SU{2) 
such that /™,„(r,z,0) = e and /„,„(r, z, 27r) = (^(2n + 3m), [/(^, 0, 0)) for all 



2 € K., r > 0. Now we define fm,n by /, 
corresponding order-parameter is 



^ir,z,ip) ^ (|(2n + 3m),C/(^,0,0)). The 



V 



73' 








(13) 



This is also the order-parameter that minimizes the vortex energy. The superfluid 
velocity and angular momentum are v = ^^{^'^ + if ) ^'^^ ^ = Nh{m + 3^)6^. 



5. Effects of external magnetic field and magnetization 

So far we have assumed that there is no external magnetic field present. Now we 
consider the situation where the condensate is placed in a magnetic field directed 
parallel to the z-axis. To the second order in the strength B of the magnetic field the 
energy from the field can be written as (23 

EbW = / d'rn{r)[jB{r){F,) + eB'{r){F^)], (14) 
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where 7 and e are constants. For ^^Rb and ■^^Na these are 7 = ±^ and e = ^j^E^f' 
Here fiB is the Bohr magneton, AEhf is the hyperfine spUtting between F — 2 
and F — I states and upper (lower) sign refers to F = 2 (F = 1). In the 
presence of a magnetic field the conservation of magnetization has to be taken into 
attention. If the magnetic field is parallel to the z-axis the magnetization is defined by 
M = J d'^rn{r){Fz) = J d^r n{r)J2f=-F J = Magnetization can have 

any value between —NF and NF, where N is the particle number. In the absence 
of an external magnetic field the magnetization is a meaningless quantity because of 
the lack of a well-defined quantization axis of the hyperfine spin. In the presence of a 
magnetic field the magnetization is a well-defined and often conserved quantity since 
the collisions that do not conserve magnetization occur usually in a timescale much 
longer than a typical lifetime of the condensate. The conservation of magnetization has 
been experimentally verified jBllHl- Thus the energy of the system should be minimized 
under the assumption of fixed M. If we consider weak values of the magnetic field it 
is enough to include only the first term in equation ifTiji in the energy. If the magnetic 
field is spatially constant this term is just magnetization multiplied by a constant. 
Thus the ground state is determined by the value of magnetization, and does not 
depend on the strength of the magnetic field, as long as the quadratic term in the 
magnetic field is neghgible. 

If the external field is absent the minimum of energy is obtained (if F = 1) 
when either |(F)| = or |(r)| = 1. If external field is present the energy has to 
be minimized under the assumption of conserved magnetization. For a ferromagnetic 
F = 1 condensate any magnetization can be produced while |(F)| = 1. Thus in 
this case the ground state is obtained if one finds an order-parameter which produces 
the given magnetization and fulfills the condition |(F)| = 1. On the other hand, 
if 1(F) I — the only possible value for magnetization is zero. If the magnetization 
is non-zero and the system is antiferromagnetic, the condition |(F)| = cannot be 
satisfied everywhere. In this case the calculation of the ground state is more difficult. 
The situation is similar in the cyclic and polar (antiferromagnetic) phases of F = 2 
system. 

Next we assume that only the term linear in the magnetic field is included 
in the energy and show that the previously calculated vortex of a ferromagnetic 
F = 1 condensate is possible for an arbitrary value of the magnetization. Direct 
calculation shows that for the vortex of equation © the magnetization is N cos^ (3. 
Any magnetization between and N can be obtained as (3 is varied between ^ and 0. 
The vortex with minimum energy is obtained when magnetization vanishes. Because 
equation Q is obtained using |F — l^mp — 1) as the reference order-parameter, 
it is not surprising that positive values of magnetization are favoured by it. A 
vortex configuration which has negative magnetization can be obtained by choosing 
|F = l,mF = — 1) as the reference spinor and calculating a representative of a vortex 
as before. The magnetization turns out to be equal to —N cos"^ f3. The vortices 
corresponding to different choices of the reference order-parameter have different 
magnetizations, but if the magnetization is allowed to vary, these vortices can be 
converted continuously into one another. Thus they are similar from topological point 
of view. 

Next we see that setting magnetization fixed makes new kind of defects stable. If 
we choose (3 — ^ in equation I^J , multiply the spinor by e*^ and redefine a a + ^, 
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This is an alternative parametrization for a ferromagnetic spinor. From the form of 
the order-parameter one might assume that a vortex is obtained if 6* = mip and a — nip 
with m, n integers. Like before, we assume that the vortices are straight and located 
on the 2:-axis. We also require that n = n{r,z),T = T(r, z). The latter assumptions 
mean that the particle density of each spin component is cylindrically symmetric||. 
The magnetization is M — J d'^r n{r) cosT{r, z), which shows that ltT5|l can produce 
any magnetization if r is chosen properly. Topologically the decay of a vortex in one 
component can now only be achieved by converting all atoms from that spin state 
into other states. From the above spinor one sees that always two components vanish 
simultaneously. When this happens magnetization is ±iV. Thus if the magnetization 
is fixed and different from ±7V, one cannot make any spin component vanish while 
keeping the magnetization fixed. If the magnetization is different from ±7V, the order- 
parameter space is 5^ x S*^ , since the phase of two components can be chosen freely. 
The dynamical stability and dynamics of these type of vortices has been studied in 
P6, 37, 38, 39, 3^. However, their topological stability following from the conservation 
of magnetization has not been pointed out before. If the order-parameter is of the 
form ltT5|l . vortices appear in each spin component separately. This means that the 
total density does not have to vanish, unless there is a vortex in every component. 
If there is one spin component without vortex the system is a coreless vortex. For 
example by choosing m = 0,n — 1 or m = l,n = —1 one gets an expression for a 
coreless vortex. 

If the magnetization is not a conserved quantity vortices with even m can be 
converted into a uniform configuration and those with odd m are equivalent with the 
winding number one vortex of equation Q. This is because in zero field r can be 
converted into a map for which r(r) = holds for all r. Then the spinor becomes 
gimip^p = l^mp = 1). This in turn is equivalent with a vortex with winding number 
zero or one for m even or odd, respectively; see 0]. 

Above we have assumed that the particle densities of different spin components 
are symmetric. It is possible that a deviation from this allows a continuous decay of 
vortex. This is also indicated by numerical studies |39j . 

If the magnetic field is strong, the term quadratic in the magnetic field has to be 
included in the energy. Also in this case the magnetization is conserved, and vortices 
which are possible when only the Hnear term is included remain topologically stable. 

6. Creation of defects 

Next we propose a method to create vortices in a ferromagnetic condensate. These 
vortices are topologiccally stable also in the absence of an external magnetic field. 
In 01] a way to create a vortex exploiting the spin degree of freedom has been 
studied. In this method a condensate in a low-field seeking state is prepared in an loffe- 
Pritchard trap. Initially the magnetic field in the z-direction is assumed to be much 
stronger than the magnetic field B± in the xy-plane. The z-component of the field 

II These kinds of vortices are called axisymmetric in 1371 1391 



we get 




(15) 



Explicit expressions for the topological defects of spinor Bose-Einstein condensates 14 

Bz is then reversed slowly, while keeping Bi_ fixed. In this way a vortex with winding 
number 2F is created, F being the hyperfine spin of the condensate atoms. The 
feasibility of this method has been experimentally verified using a ^^Na condensate 
prepared in the low-field seeking states \F = l^mp = —1) or \F = 2,mF = +2) 
|4H) . Since a vortex with a winding number larger than one is energetically unstable 
against decay into winding number one vortices, it is presumed that vortices created 
this way will split. In the case of F = 1 ^^Na condensate this has been seen to 
occur 03], and it is expected to happen also in the F = 2 case. In addition to 
creating defects which are stable in the presence of magnetic field, we propose how a 
modification of this method can be used to create 2F vortices with winding number 
one in a ferromagnetic condensate with hyperfine spin F . These vortices are stable 
in the absence of an external magnetic field. To create these vortices, in addition to 
an loffe-Prichard trap, an optical trap is needed. First a vortex with winding number 
2F is created in the previously described way. Then one waits until the vortex decays 
into winding number one vortices. Then is turned off, and the optical trap is 
turned on simultaneously. This does not change the spin state because it is assumed 
that Bz >• -B±. The remaining field in the z-direction can then be reduced to zero, 
and vortices are allowed to evolve in the optical trap. If the magnetic trap is turned 
off before the vortex has split, the vortex can in principle continuously convert into 
a non-vortex state. This is possible, for example, if the vortex reverses the rotation 
which was used to create it. This is prevented by letting the vortex split before turning 
off the loffe-Pritchard trap. In experiments there are usually stray ac magnetic fields, 
so the creation of zero-field defects requires effective magnetic shielding, which can be 
troublesome. Achieving this would, however, be rewarding, as it would allow one to 
change the topological stability of vortices as a function of the magnetic field. 

The creation of monopoles is not as straightforward as that of vortices, but a 
method for this has been proposed in |45j . 

7. Conclusions 

In this paper we have first derived a systematic way to create explicit expressions for 
vortices and monopoles. This method requires the calculation of the first and second 
homotopy groups of the order-parameter space G/H. This can be achieved by using 
equations obtained from the exact sequence of relative homotopy groups. After this 
the expressions for defects can be constructed by finding suitable mappings from the 
physical space into the group G. 

We have created examples of vortices and monopoles in spinor Bose-Einstein 
condensates using this method. Especially the defects in zero external magnetic field 
have been discussed. We have presented examples of vortices in ferromagnetic F = 1 
and F = 2 condensates and vortices and monopoles in antiferromagnetic F — \ 
condensate. We also pointed out that the order-parameter space of the cyclic phase 
of F = 2 condensate consists of two disconnected sets. The properties of one of the 
sets have been studied previously in 0. Here we calculated the topological defects of 
the other set and showed that vortices classified by Z x Zg are topologically stable, 
whereas monopoles are not possible. Also the superfiuid velocity induced by the defects 
is examined. It has been shown that in a ferromagnetic condensate with hyperfine spin 
F the presence of a vortex with winding number F does not have to induce non-zero 
superfiuid velocity or orbital angular momentum. 

We have also studied the effect of a magnetic field, concentrating on a 
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ferromagnetic F = 1 condensate. It has been shown that a vortex which is 
topologically stable in the absence of a magnetic field is also possible if the 
magnetization has a fixed value. This means that it is also possible if there is a 
homogeneous magnetic field present. In addition to this, we have found out that 
the conservation of magnetization may stabilize vortices which arc not topologically 
stable if the magnetization can vary freely. Thus there can be a transition from a 
vortex state to a non-vortex state as the magnetic field strength is lowered to zero. 
Finally a method to create vortices which are topologically stable in the absence of a 
magnetic field has been suggested. 
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Appendix A. Spin matrices and rotation operators 



The spin matrices for F = 1 are 

10 

1 1 
1 

and those of = 2 are 
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(A.2) 



\ -2 / 

Here we have set K = \. Depending on the system studied, the elements of SU(2) 
have been written either in the form 



cos 4e-*("+'^)/2 



sinfe-'(^- 



or in the form 
y(r,a,/?) = e-'*"- = 



cos ^ — z sin ^ cos (3 



.)/2 



-le 



-sinfe^(^-°)/2 
cos |e*("+T)/2 



' sin ^ sin /? 



(A.3) 



(A.4) 



— ie"* sin ^ sin /3 cos ^ + i sin ^ cos (3 

where n = (cos a sin /3, sin q sin /3, cos /3) and a = {<Jx,o'y,az) is a vector 
formed from Pauli matrices. The F-matrix has been used when discussing the 



Explicit expressions for the topological defects of spinor Bose-Einstein condensates 16 



ferromagnetic condensates, otherwise the f/-matrix has been used. 2F + 1 - 
dimensional irreducible representations of these matrices are given by the maps 
C/(a,/3,7) ^ [/(^)(a,/?,7) and V{T,a,l3) ^ V'^^^t, a, P). Here U^^\a,^,j) = 
exp{—iaFz) exp{—i(3Fy) exp{—i^Fz),V^-^\T,a, (3) = cxp(— irn- F) and F is the spin 
operator of spin F system. In this paper we need explicit expressions for U^^^ and 
V(^) for F = 1, 2. For F = 1 these are 



C/W(a,/3,7) 



e-i("+7) cos2 § -e 
ei(«-7) sin^ I 



v2 

cos/? 

■^sin/J e*("+T)cos2f 



— e*'''-4= sin/3 
^ ,2 ^ 



(A.5) 



and 



yW(T,a,/3) 



(cos f — i cos/3 sin 



V 



-\/2e'" sin/3 sin | (i cos ^ + cos f3 sin |) 
■ sin^ p 



-e^*" sin^ 



g- ia gjj^ p (-1+COS t) cos sin r 

cos2 § + COS (2/3) sin^ § 

_gia gin /3 (~ ^+'^'"" T jgji, :J-i sin r 



— e 



I a (-l+cosT) cos/3+isinT 

smp ^ 



(cos + i cos (3 sin 
For F = 2 the spin rotation matrix is 
t/(2)(a,/3,7) = 

/ e-2*(«+T) cos^ f _e-i(2«+7) sin /3 cos^ | 
g-i(a+27) sin /3 cos^ I e-'(°+T) i (cos /3 + cos 2/3) 



-2ia V6 



^sin^/3 



g-2j7x^sin2^ 



'^sin 2/3 



oi(a-27) 

e' 



4 

-e-*"^ sin 2/3 
i(l + 3cos2/3) 



e'^"-' sin /3 sin^ f e'^"-T^ i(cos /3 - cos 2/3) 



2i(a-7) sin^ I e'(2«-T) sin /3 sin^ f 



^ sin 2/3 



2*«^sin2/3 



(A.6) 



(A.7) 



-e-'(2"-7)sin/3sin2 | 

,-i(a-7) 1. 



-2i(a-7) „• 4 13 
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(cos /3 - cos 2/3) -e-^("-27) sin /3 sin^ | 



_f>«7 



87 \/6 



sin 2/3 



„2i7 V6 



sin^/3 



gi(a+7) 1 (cos /3 + cos 2/3) -e*(«+27) sin /3 cos2 | 



' sin /3 cos f 



e2i(a+7) cos4 | 
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